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ABSTRACT 


The  formulation  of  the  Griineisen  coefficient  based  on  the  velocity 
doubling  approximation  is  used  to  define  a  normalized  Griineisen  coefficient. 
A  new  integral  formulation  for  the  free  surface  velocity  is  then  written  in 
terms  of  this  normalized  coefficient.  On  the  assumption  that  the  specific 
energy  of  the  solid  at  0  "K  is  a  known  function  of  the  specific  volume  and 
that  the  bulk  sound  speed  in  the  uncompressed  state  is  a  known  quantity, 
the  Hugoniot  of  the  solid  is  chosen  to  be  that  curve,  among  a  family  of 
curves  lying  on  a  Mie-Grtineisen  constraint  surface,  which  maximizes 
the  free  surface  velocity.  A  differential  equation  for  the  resulting  Hugoniot 
is  determined  and  its  solution  is  approximated  by  the  first  three  terms  of 
a  series  expansion.  This  expansion  furnishes  a  quadratic  expression  for 
the  shock  velocity  in  terms  of  the  particle  velocity  all  of  whose  coefficients 
are  given  by  formulas  involving  physically  meaningful  quantities.  Calcu¬ 
lations  have  been  made  in  the  case  of  aluminum  and  have  been  found  to 
agree  with  the  experimental  data  out  to  340  kb  very  closely.  A  preliminary 
check  for  sodium  metal  is  also  given. 
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THE  HUGONIOT  OF  A  SOLID  DETERMINED  BY  MEANS 
OF  A  VARIATIONAL  PRINCIPLE 

I.  INTRODUCTION 

A  much  used  form  of  an  equation  of  state  for  solids  is  the  so-called 
Mie-Gr'dneisen  relation, 

p-  Pc(v)  =-^(e  -  Ec(»)j  . 

where  p^  and  E^  are  the  pressure  and  specific  internal  energy  at  0“K,  and 

T,  the  Gruneisen  coefficient,  is  a  function  of  volume  only.  One  way  of 

implementing  this  relation  for  arbitrary  (p,  v,  E)  states  involves  taking  p, 

E  to  be  known  Hugoniot  states  p  ,  E  ,  noting  that  p  =  -E1  (v),  and  taking 

for  F  one  of  the  several  formulas  available  from  the  theory  of  lattice 

dynamics  (Slater^,  Dugdale-MacDonald^,  Zubarev- Vashchenko^), 

which  express  T  in  terms  of  E£  and  its  derivatives.  The  resulting  nonlinear 

differential  equation  is  then  solved  for  E  .  With  p  ,  E  and  T  then  known 

c  c  c 

functions  of  v,  the  Mie-Gr uneisen  relation  can  be  used  to  describe  off- 
Hugoniot  states. 

The  shortcomings  and  inconsistencies  of  these  lattice-dynamical 

formulas  for  the  Gruneisen  coefficient,  and  an  account  of  what  this  implies 

in  geophysical  calculations  relating  to  the  internal  structure  of  the  earth, 

(4) 

have  been  given  by  Shapiro  and  Knopoff  .  Considerable  interest  there¬ 
fore  attends  the  possibility  of  determining  the  Gruneisen  coefficient  by  an 
independent  approach.  In  the  work  reported  here,  based  on  a  variational 
principle,  an  equation  is  derived  which  relates  the  cold  energy  function  E^ 
and  the  Hugoniot,  so  that  knowledge  of  one  yields  the  other  and  hence  the 

Gruneisen  coefficient.  This  result  is  illustrated  in  the  case  of  aluminum, 

(5) 

where  an  analytic  formulation  for  E^  due  to  McKenna  and  Pastine  is 
used  to  derive  the  Hugoniot  function.  The  latter  function  is  found  to  agree 
with  available  experimental  data  very  closely  up  to  340  kb.  An  interesting 

1 


V 


4 


,*■  rmHt&nxmil  Wi^MfUlNW^B 


by-product  in  this  regard  is  a  quadratic  expansion  in  the  particle  velocity, 
which  approximates  the  shock  velocity;  the  coefficients  of  this  expansion 
are  simple  formulas  involving  physically  meaningful  quantities. 
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II.  ANALYSIS 

As  was  shown  recently  by  Heyda^  the  temperature-independent 
Griineisen  coefficient  of  a  solid,  T (u),  may  be  expressed  in  the  form 


r(u)=r*(u)G  , 


(1) 


where  F*(u)  is  the  form  which  F  takes  on  the  basis  of  the  approximate 
doubling  of  the  particle  velocity  u  behind  a  planar  shock  wave  moving 
through  the  solid  with  speed  D  and  arriving  at  a  free  surface  parallel  to 
the  shock  front.  The  function  r*(u)  was  found  to  be 


r*(U)  _  (D-U)(D^D')D' 

D2 


(2) 


where  the  prime  indicates  differentiation  with  respect  to  u.  The  function 
G,  which  can  be  regarded  as  a  normalized  Griineisen  coefficient, 


G  =  T(u)  /  T*  (u)  , 


(3) 


was  shown  in  turn  to  be  given  by 


(4) 


where  u  =  u,  -  u  and  u.  is  the  free  surface  velocity. 

r  fs  fs  7 


Solving  Eq.  (4)  for  u'  and  integrating  with  respect  to  u,  we  may 

r 


write 


-  e 


2  1 


e2G 


du  , 


(5) 


where 


e  =  uD  /  D 


(6) 
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It  is  convenient  to  change  tiiv  .^/dependent  variable  from  a  to  c  ai.,1 

to  regard  u,  G,  u  and  D  as  functior  •  of  <r..  .If  c  shou'd  be  know:  from  „.n 
r 

independent  source  as  a  function  of  u,  t*. e  relation  may  be  looker,  up 
as  a  differential  equation  for  the  determii  *.tio  :  of  the  Tfugoiiiot  function 
D(u).  Making  this  change  we  can  rewrite  Eq,  {5>  '  tkr  form 


r£  . 

u  =  /  u 

r  Jn 


1  -  e  G(e) 


where  now  the  prime  indicates  differentiation  with  respect  to 

For  e  arbitrary  but  fixed  it  is  natural  to  regard  s  u  u]  as  ?. 
functional  defined  on  the  set  of  functions  G  and  u  and  to  consider  e-<tr  e'nizing 
u  [G,  u]  over  a  family  of  curves  G  =  G(e),  u  =  u(e)  lying  on  some  cc;/~ 

r 

straint  surface  g  a  g(e,  G,  u)  =  0. 


To  obtain  a  suitable  function  g  we  refer  to  the  Mie-Gruneisen  form 
of  the  equation  of  state  of  a  solid,  namely, 


P  -  Pc(v) 


=  [E-  Ec<v>]  • 


where  v  is  specific  volume  and  p,  E  denote  pressure  and  specific  energy, 
respectively;  the  subscript  c  indicates  pressure  and  specific  energy  at 
0°K.  Taking  p,  E  along  the  Hugoniot  of  the  solid  centered  at  the  normal 


conditions 


pH  =  °'  EH  =  Ec  (vo>  +  Eo  ’ 


where  v  is  the  specific  volume  of  the  solid  at  T  =  T  ,  p  =  0,  and  E  is 
o  1  o  r  o 

the  specific  thermal  energy  of  the  solid  at  these  same  conditions,  (in  our 
calculations  of  the  aluminum  hydrostat  we  shall  take  Tq  =  300  °K),  we 
may  write 


—  [  E  -  E  1  -  p  +p=0  . 

v  L  c  H  J  rc  rH 
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1 


Since 


and 


o  -  -  a  „  (v) 
c  c 


(11) 


eh  =  eh  i_r|v0-v)  ■ 


(12) 


can,  upon  making  the  definition 


£v  ,v)  (v)  .  E  (v  j  .  E  , 

C  !  c  I  c  o  o 


(13) 


rewrite  Eqj  (10)  as 


-v,1 

v  L  c  2  o  J 


(Y0- V)J +ec(v)+Ph  =  0  • 


(14) 


^rom  the  Hugoniot  relation  (for  a  derivation  of  the  Pankine-Hugoniot 
relations),  see  Ref.  (7)) 


v  =  v 


(>-s) 


1 


(15) 


we  obtain 


dv 

du 


vo(l-€) 


(16) 


so  that 


dE 

_ c 

dv 


dE 


v  ( 1  -  e)  du 
o 


or,  in  a  less  accurate  but  more  convenient  notation, 


E’  (v)  =  - - £ - -  E'  (u)  . 

c  v  (1  -  e)  c 


(17) 


Recalling  that 


5 


1 


•  \ 


PH  =  PoDu'  PH(vo  -  v>  =  “  •  Po=  I/vo  ' 


we  find,  upon  substituting  from  Eqs.  (15),  (17),  (18)  into  Eq.  (14),  that 


d^t  1E=  "2“  J -ttv  +  “  =  °  • 


From  Eqs.  (1)  and  (2),  however,  it  follows  that 


e(l+e)G(e)/u  , 


which,  substituted  into  Eq.  (19),  yields  the  desired  form  of  the  constraint 
surface  g  =  0,  namely. 


g(e,  G,  u)  =  e(l  +  e)G(e)-rivj,  e)  =  0  , 


where  the  function  H  is  defined  by 


In  the  space  of  the  variables  e,  G,  u  the  equation  (21)  constitutes  a 

surface  passing  through  the  curve  G  =  G  (e),  u  =  u(e).  For  each  choice  of 

the  function  G,  there  is  determined  a  corresponding  function  u(e). 

Assume  then  that  the  functional  u  [G,  u]  has  an  extremum  for  the  curve 

r 

l  j 

G  =  G  (e),  >  u  =  u(c)  .  (23) 


Then  since  g_  and  g  do  not  vanish  simultaneously  at  any  point  of  the 
G  u 

surface  (21),  there  exists  (see  Theorem  2,  p.  46  of  reference  (8))  a 
function  X(e)  such  that  the  curve  (23)  is  an  extremal  of  the  function 


f  [F  +  X(e)  g  ]  de  , 


\ 

i 


where  F  is  defined  oy 


MPWWNSWWwijp  'few 


F  =  a  (e) 


1  -  e 

1  -  e2G  (e) 


the  necessary  conditions  for  which  are  given  by 


FG,XgG  =  0-  X6u  =  d7Fa'  • 


Eliminating  X  in  Eqs.  (25)  and,  for  convenience,  letting 


2 

1  -  e  G 


we  obtain  the  single  equation 


u'  P  H  =  €  ( 1  +  e)  P  +  g'p^ 
G  u  e  G 


However,  from  Eq.  (21)  we  have 


e  ( 1  +  e)  G*  +  ( 1  +  2e)  G  =  H  u'  +  H  , 

u  e 


so  that  u'  and  g'  can  be  eliminated  simultaneously  from  Eqs.  (27),  (28) 
and  no  differential  equation  need  be  solved  to  obtain  G(e)  and  u(e).  The 


result  is 


(1  +  2e)G 


-J 


PG  =  e  (1  +  e)Pe  . 


When  Eq.  (26)  is  taken  into  account  in  Eq.  (29)  and  the  resulting  differen¬ 
tiations  carried  out,  we  find  the  result 


2  -  (1  -  e)H 


G  (e)  = 


1  -  e  +  2e 
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From  Eq.  (22)  we  note  that 


u 


(1  -  e)  H  =  H  +  •  ,  , 

e  *  1  2 

E  -  -  u 
c  2 


=  e  (1  +  c)  G  + 


_*  1  2 

E  -  -  u 
c  2 


so  that  Eq.  (30)  assumes  the  final  form 


G(e) 


2 


1  +  3e 


2 


(31) 


where  the  function  u  =  u(e)  is  obtained  by  eliminating  G(e)  between  the 
Eqs.  (21),  (31),  the  function  E*  being  assumed  to  be  known.  Actually,  it 
is  easier  to  obtain  the  function  e  =  e  (u).  We  find  that  it  satisfies  the 
equation 


3 

e 


2 

+  u 


+  3c 


"u  E  (u)  -  u  +e  3u^  -  2E*  +u  E  (u)-u  =0.  (32) 

c  [  c  l  c 


From  the  relation  (17)  we  may  readily  show  that 

uE^(u)  =  (1  -  e)(x  -  xq)E^(x)  ,  (33) 


where  x  =  v/v  and  v  is  the  specific  volume  at  T  -  0  °K  and  p  =  0. 
oK  o  K 

With  the  relation  (33)  we  may  then  recast  the  cubic  equation  (32) 
in  the  form 


,  2  * ...  .2 
(  u  -  Ec)  (1  -  e) 


=  (1  +  3c  )  (x 


x  )  E  (x)  -  E* 
o  c  c 


(34) 


where  a  factor  (1  -  c)  has  been  divided  out.  Equation  (34)  can  now  be 
solved  for  e,  and  hence  for  d',  to  get 


8 


I 

I 


0 


* 


where 


c 


(35) 


Since 


c 


P 


c 


and 


x  =  x  ( 1  -  ~) 
o'  D' 


x 

o 


V 

o 


(36) 


we  may  rewrite  Q  as 

v  up 
o  c 

Q  =  _ — 

U2  -  Ec  4  E„ 

o 

upon  using  the  relations  (13)  and  (9). 

We  note  that  Q  -*  1  as  u  -» 0.  Hence  the  minus  sign  in  Eq.  (35)  must 
be  chosen  in  order  that  d'  remain  finite.  The  indeterminacy  in  the  right 
member  of  Eq.  (35)  at  u  =  0  may  be  resolved  by  rewriting  it  in  the  form 


E  + 

c 


EHr 


D 

u 


1  -  Q 

1  +  V  4 Q  -  3Q2 


and  noting  that 


1  -  Q 
u 


u  - 


V  p 

o  c 


Ec  +  EH 
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Since  E£  and  p^  are  generally  known  as  functions  of  x,  and  x  in  turn 
may  be  expressed  in  terms  of  u  and  D  by  means  of  relation  (36),  we  see 
that  the  Hugoniot  function  D  =  D(u)  may  be  obtained  by  solving  the  first 
order  nonlinear  differential  equation 

D'  =  /( u,  D)  ,  (37] 


where 


/  (u,  D)  = 


R(u,  D) 


1  +  V4Q(u,  D)  -  3Q2(u,  D)' 


V  u  p 
o  c 


-  E  +  Eu 
c  «0 


“  -Ec  +  EH 


uD  -  v  p 
o  c 


□  -ec  +  eh 


subject  to  the  initial  condition 


u  =  0  ,  D  =  D (0)  =  c 

o 


Since  a  closed  form  solution  of  Eq.  (37)  appears  unlikely,  the 
solution  D(u)  must  be  obtained  either  through  numerical  integration  or 
expansion  in  series.  If  we  denote  the  expansion  for  D  through  the 
quadratic  term  by 


D  =  c  +su  +  bu 
o 


we  find  by  putting  u  =  0  in  Eq.  (38)  that  s  is  given  by  the  simple  formula 

v  p  (v  ) 
o  c  o 


or  equivalently, 


10 


8 


X 

o 


E  (X ) 
c  o 

2E 


(43) 


Continuing,  we  find  upon  differentiating  (37)  with  respect  to  u  and  then 
putting  u  =  0  that 


b  = 


5L 1121  = 


Z  it  .2 

x  E  (x  )  -  6s  E 
0  c  o  o 

4E  c 
o  o 


(44) 


It  might  be  noted  that  the  value  of  a  given  by  formula  (43)  will  normally 

differ  from  tabulated  values  of  the  slope  of  the  "linear"  shock  velocity- 

particle  velocity  relation  since  the  latter  represents  an  empirical  fit  to 

experimental  data  over  a  range  of  particle  velocities  well  beyond  u  =  0. 

TheHugoniot  formulation  (38)  will  in  addition  reflect  passage  through  a 

minimum  point  on  the  curve  E  =  E  (v)  at  v  =  v  as  v  decreases  from 
r  c  c  oK 

v  ,  i.  e.  as  u  increases  from  u  =  0. 
o 


Additional  coefficients  beyond  D  (0 )/2  may  be  calculated  by  the 
same  procedure.  However  existing  analytical  formulations  for  E  (x) 
do  not  offer  sufficient  accuracy  in  the  values  of  Ec'  '  (x^)  beyond  n  =  2. 
Anticipating  that  such  formulations  may  become  more  accurate  in  the 
future  we  have  calculated  ore  additional  coefficient: 


h  p 


D"'(0)  x  3  e'"  (x  )  -  24sE  c  d"  (0)  -  8E  s3  +  2s c  2 
_  _  o  c  o _ o  o _ o _ o 

6  1 2  E  c  2 

0  o 


(45) 


The  Griineisen  coefficient  T  (u)  is  now  readily  available  from 
Eqs.  (2),  (3),  and  (31).  We  find 

/ 

2 


r(u)  = 


2  [  D*  (u)  -  e]  (l+€) 


1  +  3e 


:  +u2  -  Te  -  E  (v  )1 
_o _ l  c  C  O  J 

7  u2-[e  -  E  (V  )1 
2  L  c  c  o  J 


(46) 


E  + 
o 


11 


or,  equivalently, 


Hu)  =  (D  -u) 


e'  (u) 

C 


1  -  G 


(47) 


The  initial  value  of  T  turns  out  to  be  T  (0)  =  2s,  whence  from  Eq.  (43)  we 
obtain 

x  E*  (x  )  /  \ 

r  (0)  =  2s  =  — ^C— -  ,  f  G  (0)  =  2  j  ,  (48) 

a  result  which  follows  directly  from  the  Mie-Grtineisen  relation  (14)  by 
putting  v  =  v  . 

That  the  extremal  (23),  where  G(e)  and  u(e)  are  defined  by  the 
relations  (31)  and  (32),  maximizes  the  functional  u^tG,  u]  is  evident  from 
the  fact  that  for  the  choice  G  s  1  (which  corresponds  to  taking  T  in  accord¬ 
ance  with  the  velocity  doubling  approximation)  we  obtain  u^  =  u,  whereas  a 
three  term  expansion  of  the  integral  in  Eq.  (5)  yields  the  result 


u 

r 


u  + 


■tno)-.]  u3 

6c 

o 


(49) 


which  in  view  of  the  result  (48)  shows  that  u  >  u.  Indeed,  with  c  assumed 

r  o 

to  be  known  and  with  Dsc^t  su,  the  D  function  (or,  equivalently,  the  G 
function)  which  maximizes  uf  is  that  for  which  s[T(0)  -  s  ],  considered  as 
a  function  of  s,  is  a  maximum.  This  occurs  evidently  for  T (0)  =  2p,  in 
agreement  with  relation  (48). 


The  approximate  relation  (49)  is  of  interest  in  its  own  right  as  it  gives 

good  agreement  with  tabulated  values  of  u  /u  for  many  of  the  metals  studied 

(7)  r  ' 

by  Rice,  McQueen,  and  Walsh  .  Using  the  T  and  s  values  listed  there, 

o 

we  find  the  comparisons  for  copper,  24  ST  aluminum,  silver,  and  zinc 
shown  in  Table  I. 


TABLE  I 


COMPARISON  OF  VALUES  FROM  EQ.  (49) 
WITH  THOSE  GIVEN  IN  REF.  7 


Copper 


p(kb) 

u  (cm/usec] 

|  ur/u[Eq.  (49)] 

ur/u  [Ref.  7] 

100 

0.0263 

1.001 

1.001 

300 

0. 0681 

1.004 

1.005 

500 

0. 1025 

1.010 

1.012 

24  ST  Aluminum 

100 

0.0580 

1. 002 

1.003 

300 

0. 1465 

1.013 

1.015 

500 

0.  2145 

1.028 

1.030 

Silver 

100 

0.0256 

1.  002 

1.001 

300 

0.0658 

1.011 

1.010 

500 

0.0992 

1.024 

1.022 

100 

0.0378 

Zinc 

1.003 

1.004 

300 

0.0925 

1.019 

1.022 

500 

0.  1373 

1.043 

1.042 
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III.  EXAMPLE  FOR  ALUMINUM 


For  aluminum  we  employ  the  analytical  formulation  given  by 
McKenna  and  Pastine  in  Ref.  5, 


Ec  <x|  =  VoK 


(t)[ 


2(3(1 -x1/3) 


-  2e 


P(l-x1/3) 


I 


(50) 


where 


v  .  =  0.  366  cc/gm,  v  =  0.  371  cc/gm, 
oK  o 

A  =  0.6271  Mb,  0  =  3.  772 

E  =  .  00164  Mb  cc/gm,  c  =  .  5404  cm/nsec,  T  =  300  #K. 


Substituting  into  formulas  (43),  (44),  and  (48)  we  obtain 

8=1.  171,  b  =  1.  303  jjsec/cm,  T  =  2.  342, 

o 

so  that  the  Hugoniot  function  D  is  approximated  by  the  three  term 
expansion 

D  =  0.  5404  +  1.  171  u  +  1.  303u2  .  (51) 


This  function  is  plotted  in  Figure  1  and  agrees  with  the  experimental 

data  for  24  ST  aluminum,  given  by  Walsh  and  Christian  in  Ref.  (9)  t-.d 

Walsh,  Rice,  McQueen  and  Yarger  in  Ref.  (10),  remarkably  well  up  to 

u  =  0.  165  cm/usec,  which  corresponds  to  a  pressure  of  approximately 

340  kb.  Byond  this  point  the  curve  diverges  slowly  from  the  data  taken 

from  Ref.  (10).  An  additional  term  in  Eq.  (51)  would  remedy  this  if 

derivatives  of  E  (x)  of  the  third  and  higher  orders  were  known  at  x  =  x 
c  o 

with  sufficient  accuracy.  This,  however,  puts  too  great  a  strain  on  the 
empirical  form  (50). 


14 


V 
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IV.  CONCLUSION 


The  tentative  conclusion  is  that  the  aluminum  hydrostat  is  deter¬ 
minable  from  the  variational  principle  proposed  here  and  that  the  same 
may  well  be  true  for  other  solids.  Accepted  as  a  general  principle,  the 
procedure  may  be  reversed  and  the  cold  energy  may  be  obtained  by 
solving  the  linear  differential  equation  (32)  for  with  an  accurate  experi¬ 
mentally  determined  Hugoniot  as  input.  Calculation  of  the  Griineisen 
coefficient  can  then  be  made  directly  from  the  Mie-Griineisen  relation  (8), 
thereby  avoiding  any  use  of  formulas  based  on  the  theory  of  lattice 
dynamics. 


0©ST//UJO 

o  A 110 013 A  XOOHS 


J. 

CM 

o 


go 

c 

u 

u 

C 

• 

«M 

3 

3 

m 

O 

h- 

O 

m 

• 

V 

-a 


jo 

£ 

Q 


W 

H  (J  +* 

•  «  2 

+  ■§  0 

2  » 


3 

> 


o 

m 


xtJC 


°  £ 


y  .2 

2  x> 


cS" 

>■3- 

uig 


cc 

s 


II 


y  w 

>  to 

o 

A 


U  a» 

*>  2 
O  lllH 
■H  r* 

C  U  <t 

O  2  in 
»  «  0 

X  g  ■ 

>— 1  «i 
Oi  ^  rZ 

o  <0  V 

•IS  A  <o 

4>  **  X 
£  DC  +» 

0  C  t. 
0)  -  U 

43 


« s 

■rl  ■> 


<U 

43 


•5  C 

»  »  ° 

2  g.3 

C  u  T> 
H>  !»C 
(0,0  0 
0)  u 

L  1 0 

E  di  +> 

g  c 

»-  ’2  « 


0) 

> 

u  _ 

3  a>  c 

Q  T3  C 


0) 


0) 


4J  C 


V 

u 

3 

DO 


■ 


16 


REFERENCES 


Slater,  J.  C.  ,  Introduction  to  Chemical  Physics,  Chapter  XIII, 
McGraw-Hill,  New  York  (1939). 

Dugdale,  J.  S.  and  D.  K.  C.  MacDonald,  "The  Thermal  Expansion 
of  Solids,  "  Phys.  Rev. ,  _89,  832  (1953). 

Vashchenko,  V.  Ya.  and  V.  N.  Zubarev,  "Concerning  the  Griineisen 
Constant,"  Soviet  Physics  -  Solid  State,  5, ,  No.  3,  653  (1963). 

Shapiro,  J.  N.  and  L.  Knopoff,  "Comments  on  the  Interrelationships 
between  Griineisen's  Parameter  and  Shock  and  Isothermal  Equations 
of  State, "  J.  Geophys.  Res.,  74,  No.  6,  1439  (1969). 

McKenna,  P.  and  D.  J.  Pastine,  "Volume  Dependence  of  the 
Griineisen  Parameter  for  Aluminum,  "  J.  Appl.  Phys.  _39,  6104 
(1968). 

Heyda,  J.  F.  ,  "The  Griineisen  Coefficient  Implied  by  the  Velocity 
Doubling  Rule,  "  J.  Appl.  Math,  and  Phys.  (ZAMP),  ZZ,  Fasc.  2, 

350  (1971). 

Rice,  M.  H.  ,  R.  G.  McQueen  and  J.  M.  Walsh,  "Compression  of 
Solids  by  Strong  Shock  Waves,  in  Solid  State  Physics,  Vol.  6, 
Academic  Press,  New  York  (1958). 

Gelfand,  I.  M.  and  S.  V.  Fomin,  Calculus  of  Variations,  Prentice- 
Hall,  Englewood  Cliffs,  New  Jersey  (1963). 

Walsh,  J.  M.  and  R.  H.  Christian,  "Equations  of  State  of  Metals 
from  Shock  Wave  Measurements,  "  Phys.  Rev.,  97,  1544  (1955). 

Walsh,  J.  M.  ,  M.  H.  Rice,  R.  G.  McQueen  and  F.  L.  Yarger, 
"Shock-Wave  Compression  of  Twenty-Seven  Metals.  Equations  of 
State  of  Metals,  "  Phys.  Rev.,  108,  No.  2,  196  (1957). 

Pastine,  F.  J.  ,  "Equation  of  State  for  Sodium  Metal,  "  Phys.  Rev.  , 
166,  No.  3,  703  (1968). 


/ 


APPENDIX 

We  list  here  some  additional  formulas  of  interest  which  are  derivable 
from  those  obtained  in  Section  II. 

First  of  all  we  may  readily  show  that  Eq.  (32)  may  be  converted  to 
the  form 


1  -  e  +  2e  - 


2  2 


_2i£L 


1  +  3e 


G(e) 


H 


(A  -  1) 


which  relates  the  Hugoniot  and  the  cold  pressure  in  an  interesting  way. 

Next  by  requiring  this  relation  to  hold  in  the  limit  as  e  -»  0,  we  obtain 


G'(e) 


=  0  . 
e=0 


G"  (£) 


2c 


=  -12  + 

e=0 


E  s 
o 


2 


(A -2) 


and  the  limiting  result 


Pc(V  = 


4p  c 
o  o 


12  +  G  (e) 


L 


-  2  p  E  s 
o  o 


(A-  3) 


From  Eq.  (4),  however,  we  find  that 


G  (0) 


(«)  (0)  _ 

s  ^ 


so  that  in  conjunction  with  relation  (A-2)  we  obtain  an  additional  term  for 
the  series  (49),  namely, 
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M 


/  n  2  4 

T-t-  c  D  (0)  -  s  a 
3  o 


(A-4) 


Finally,  by  differentiating  Eq.  (3),  taking  note  of  Eq.  (2),  and 


again  employing  relation  (A-2),  we  find 


r'  (v  )  =  -f-  [  c  d"  (0)  +  s  (s\-  1)  ]  . 


O  V  o 
o 


(A -5) 


Close  to  x  =  x  ,  then,  we  have 
°  i 

1  r  (X)  =  r  (x  )  +  r ' 


(x)  =  r  (x  )  +  r  (x  )  •  (x  -  x  )  ,  x  =  v  /v 

o  o  o  o  OK 


Thus  for  aluminulm, 


r (1>  =  r (1.  01366)  +  3. '714(1  -  1.01366) 
=  2.  300 


This  compares  favorably  with  McKenna  and  Pastine'  ,  who  find 

r  (1)  =  2.301.  \ 

A  final  example  is  the  evaluation  of  T  =  2s  for  metallic  sodium. 

o 

From  Eq.  (43|  we  have 

I  1 

XoE  c(XQ)  (1.0299)  (0.00314)  ,  w 

.  V  El  =  0.002784  =  1*  16  ’ 

I  ° 

which  compares  weUwith  the  value  of  1.  15  obtained  by  Pastine  in  Ref.  11, 
in  which  he  evaluated  the  role  of  anharmonic  contributions  in  determining 
T  ty  a  lattice  dynamical  approach.  In  this  calculation  we  obtained 
Pc  =  -.  00314  Mb  by  interpolating  in  Table  I  of  Ref.  (11).  The  value  of  E 
was  obtained  from  the  Debye  expression 


E  =  r  D 
o  A  i 


, 

1 V  T„  /  ’ 


I 
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where  R  =  8.  314  x  10  Mb  cm  /deg,  T  =  300  #K,  A  =  22.997  gm,  and 

o 

Dj  is  the  Debye  function  with  argument  6(XQ)  =  140  °K  obtained  from  an 
approximate  formula  given  by  Pastine^^.  The  value  of  D  j(4~)  turns 
out  to  be  0.  85578.  , 
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